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Abstract 



Some new exact solutions of Einstein's field equations have come forth within the scope of a spatially 
homogeneous and anisotropic Bianchi type-Ill space-time filled with barotropic fiuid and dark energy by 
Q , considering a variable deceleration parameter. We consider the case when the dark energy is minimally 

^ ■ coupled to the perfect fluid as well as direct interaction with it. Under the suitable condition, the anisotropic 

C/2 ' models approach to isotropic scenario. We also find that during the evolution of the universe, the equation 

^ , , of state (EoS) for dark energy ijO^'^''\ in both cases, tends to —1 (cosmological constant, lo^'^'^^ = —1), by 

displaying various patterns as time increases, which is consistent with recent observations. The cosmic jerk 
parameter in our derived models are in good agreement with the recent data of astrophysical observations 
^ , under appropriate condition. It is observed that the universe starts from an asymptotic Einstein static era 

' and reaches to the ACDM model. So from recently developed Statefinder parameters, the behavior of different 

stages of the universe has been studied. The physical and geometric properties of cosmological models are 
also discussed. 

^ ■ PACS numbers: 98.80.Es, 98.80-k, 95.36.-t-x 

Key words: Bianchi type-Ill models. Dark energy, Variable deceleration parameter, Accelerating universe 
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1 Introduction 

There is good evidence that a mysterious form of dark energy (DE) accounts for about two-third of matter and 
energy in the Universe. The direct evidence comes from distance measurements of type la supernovae (SNe la) 
as standard candles which indicate the expansion of the universe is speeding up, not slowing down [T] — [3]. In 
addition, measurements of cosmic microwave background [3] and the galaxy power spectrum also indicate the 
existence of the dark energy. These observations have reopened the quest for the cosmological constant which 
was introduced by Einstein [B] in his field equations, but later abandoned [7] and infamously cited as his greatest 
blunder [S]. The cosmological constant can be considered as new kind of "world matter" and be identified 
with the energy density of the vacuum |10j . The simplest candidate for dark energy is the energy density of the 
quantum vacuum (or cosmological constant) for which p — ~p. However, the inability of particle theorists to 
complete the energy of the quantum vacuum - contributions from well understood physics amount to 10^^ times 
critical density - casts a dark shadow on the cosmological constant [llj . In addition to this, a number of viable 
models for DE have been fabricated. These scenarios include, quintessence [Tll[l3], modified gravity [H] — [ID], 
tachyon [5T] arising in string theory [32] , quintessential inflation [53] , chaplygin gas as well as generalized chap- 
lygin gas [23] — [22]) cosmological nuclear energy [5S], equation of state (EoS) parameter [5S] — [32], braneworld 
[40[ |4T] and interacting dark energy models [42] — [48] . Therefore some form of dark energy whose fractional 
energy density is about fl^'^'^^ = 0.70 must exist in the Universe to drive this acceleration. This fact can be put 
in agreement with the theory, if one assumes that the Universe is basically filled with so-called dark energy. 
Evolution of the equation of state (EoS) of dark energy uj^'^'^'f = transfers from lu^'^'^^ > — 1 in the near 
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past (quintessence region) to u}^^'^^ < — 1 at recent stage (phantom region) 02] — [SI]- So another cosmological 
coincidence problem may be proposed: why oj^'^'^^ = 1 crossing is occurred at the present time (55]. When SNe 
resuhs are combined with five-year WMAP, it is found that -1.38 < w^''^) < 0.86 [53]-[55]. For recent review, 
the readers are advised to see the references of Padmanabhan Copeland et al. [57], Perivolaropoulos [55] . 
Jassal et al. [SS] and Miao et al. [50] . 

The simplest model of the observed universe is well represented by Friedmann-Robertson- Walker (FRW) 
models, which are both spatially homogeneous and isotropic. These models in some sense are good global ap- 
proximation of the present-day universe. But it is also believed that in the early universe the FRW model does 
not give a correct matter description. The anomalies found in the cosmic microwave background (CMB) and the 
large structure observations stimulated a growing interest in anisotropic cosmological models of universe. Ob- 
servations by the Differential Microwave Radiometers (DMR) on NASA's Cosmic Background Explorer (COBE) 
registered anisotropy in various angle scales. Ellis |61) pointed out that although the observed universe seems 
to be almost isotropic on large scales, the early and/or very late universe could be anisotropic. Furthermore, 
the interest in such models was encouraged in recent years due to the argument that going on the analysis and 
the interpretation of the WMAP [S2] — [SI] data, whether they need a Bianchi type morphology to be explained 
successfully jSS] — [70]. Thus, in recent years Bianchi universes have been gaining an increasing interest of obser- 
vational cosmology. 

In earlier, Xin |71j studied an interacting two-fluid scenario for quintom dark energy. Xin-He et al. [72j 
considered Friedmann cosmology with a generalized EoS and bulk viscosity to explain DE dominated universe. 
Recently, several authors [73] — [85] have examined and discussed the DE models in different context of use. 

Incited by above discussions, in this paper, we purport to probe the evolution of the dark energy parameter 
within the framework of a Bianchi type-Ill cosmological model filled with two fluids by considering a variable 
deceleration parameter (DP). In doing so we consider both non-interacting and interacting cases. This paper is 
organized as follows: the metric and the field equations are presented in Sec. 2. Section 3 deals with the exact 
solutions of the field equations. In Sec. 4 we obtain the non-interacting two-fluid model and its geometric and 
physical signiflcance are discussed. Section 5 deals with cosmic jerk parameter. In Sec. 6 interacting two-fluid 
scenario is presented. Concluding remarks are given in Sec. 7. 

2 The Metric and Field Equations 

We consider the Bianchi type-Ill metric as 

ds^ = ~dt^ + A^{t)dx^ + B^it)e-^°"'dy^ + C\t)dz^, (1) 
where A{t), B{t) and C{t) are functions of time only. 

We define the following physical and geometric parameters to be used in formulating the law and further in 
solving the Einstein's field equations for the metric ([1}. 

The average scale factor a of Bianchi type-Ill model ([!]) is defined as 

a^{ABC)i. (2) 

A volume scale factor V is given by 

V = a^ = ABC. (3) 
We define the generalized mean Bubble's parameter H as 

H=^{H,+Hy + H,), (4) 

where Hx = ^, Hy = g and = ^ are the directional Bubble's parameters in the directions of x, y and z 
respectively. A dot stands for differentiation with respect to cosmic time t. 
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From Eqs. we obtain 
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The physical quantities of observational interest in cosmology i.e. the expansion scalar 0, the average anisotropy 
parameter Am and the shear scalar are defined as 
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(7) 
(8) 



where A7?i = Hi — H(i = x, y, z) represents the directional Hubble parameter in the direction of x, y, z respec- 
tively. Am = corresponds to isotropic expansion. 



The Einstein's field equations ( in gravitational units SttG = c = 1) read as 



-T, 



(m)i 



{de)i 



(9) 



where T^j-™-'* and T^'^'^^^ are the energy momentum tensors of perfect fluid and DE, respectively. These are given 

by 



(m)z 



diaghp('"\p('"\p('"\p(™)] 



= diag[-f,w 



(10) 



and 



T^'^")' = diaghp^'^'^) r,('i«) r^ide)^ 



= diag[-l,a;(''^),Lj(''^\tj('''=V^'^'\ (H) 
where p^™) and are, respectively the energy density and pressure of the perfect fluid component or ordinary 
baryonic matter while / p{m) is its EoS parameter. Similarly, p^'^'^^ and p^'*'^) are, respectively the 

energy density and pressure of the DE component while I p^'^'^^ is the corresponding EoS parameter. 

We assume the four velocity vector = (1,0,0,0) satisfying = —1. 

In a CO- moving coordinate system (m* ~ 6'q), Einstein's field equations ^ with (flUl) and (|lip for Bianchi 
type-III metric ([T]) subsequently lead to the following system of equations: 
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= pi"') + p(de) ^ 
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The law of energy-conservation equation (T.y — 0) yields 

+ 3(1 + uji"''>)pi"''>H + pi'^") + 3(1 -I- uji'^''))pi'^'')H = 0. 
The Raychaudhuri equation is found to be 



= -(l + 3^(^'=))p(^'=)-i( 



2a\ 



(12) 

(13) 
(14) 
(15) 
(16) 

(17) 
(18) 



3 



3 Solutions of the Field Equations 



The field equations (|12l) - (jl6l) are a system of five linearly independent equations with seven unknown parameters 
A, B, C, p^"), p^'^''\ p^'^''\ oj^<^''\ Two additional constraints relating these parameters are required to obtain 
explicit solutions of the system. 

B = (19) 



Eq. p6]) . obviously leads to 



where is an integrating constant. 

Firstly, we assume that the scalar expansion 6 in the model is proportional to the shear scalar. This condition 
leads to 

A = C", (20) 

where n is a constant. The motive behind assuming this condition is explained with reference to Thorne j86| . 
the observations of the velocity-red-shift relation for extragalactic sources suggest that Hubble expansion of the 
universe is isotropic today within « 30 per cent |871 188j . To put more precisely, red-shift studies place the limit 
^ < 0.3 on the ratio of shear a to Hubble constant H in the neighbourhood of our Galaxy today. Collins et 
al. |89| have pointed out that for spatially homogeneous metric, the normal congruence to the homogeneous 
expansion satisfies that the condition ^ is constant. 

We define the deceleration parameter q as 

d / a\ ^ aa ( H + \ , , ^ 

= b{t) say, (21) 



a \ c I \ JJ2 



where a is the average scale factor of the universe defined by Eq. ([2|) and the dots indicate derivatives by proper 
time. The expansion of the universe is said to be "accelerating" if d is positive (recent measurements suggest 
it is), and in this case the DP will be negative. The minus sign and the name "deceleration parameter" are 
historical; at the time of definition q was thought to be positive, now it is believed to be negative. Recent 
observations [1] — [3] have suggested that the rate of expansion of the universe is currently accelerating, perhaps 
due to dark energy. This yields negative values of the DP. 

The motivation to choose such time dependent DP is behind the fact that the universe is accelerated expan- 
sion at present as observed in recent observations of Type la supernova [I] — [3] and CMB anisotropies [SSllMllOTj 
and decelerated expansion in the past. Also, the transition redshift from deceleration expansion to accelerated 
expansion is about 0.5. Now for a Universe which was decelerating in past and accelerating at the present time, 
the DP must show signature fiipping ^92| — ^94j- So, there is no scope for a constant DP at the present epoch. 
So, in general, the DP is not a constant but time variable. 



The equation (PT|) may be rewritten as 



a" 



= 0. (22) 
a 

In order to solve the Eq. (P^ . we assume h — b{a). It is important to note here that one can assume 
b — b(t) — b{a{t)), as a is also a time dependent function. It can be done only if there is a one to one 
correspondences between t and a. But this is only possible when one avoid singularity like big bang or big rip 
because both t and a are increasing function. 

The general solution of Eq. ([^^ with assumption b = b(a), is given by 

e/^'^'^da^t + fc, (23) 
where k is an integrating constant. 
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One can not solve Eq. (|23p in general as b is variable. So, in order to solve the problem completely, we have 
to choose / ^da in such a manner that Eq. be integrable without any loss of generality. Hence we consider 

-da = lni(a). (24) 
a 

which does not affect the nature of generality of solution. Hence from Eqs. (|23p and ([24| . we obtain 



J L{a)da = t + k. (25) 



Of course the choice of i(a), in Eq. (|25p . is quite arbitrary but, since we are looking for physically viable models 
of the universe consistent with observations, we consider 

L{a) = ^=L=, (26) 

where a is an arbitrary constant. In this case, on integrating, Eq. (I25|) gives the exact solution 

a(t) sinh(/3r), (27) 

where T — t + to, to and /3 being constants of integration. We also note that t + and T — oo respectively 
correspond to the proper time t = —to and t = oo. The relation ([ST)) is recently used by Pradhan et al. [38] and 
Amirhashchi et al. [47] in studying dark energy models in Bianchi type-VIo and FRW space-time respectively. 

Now, by using ([S]), ((THj), (HOj) and we can find the metric components as 

A = iismh^{(3T), (28) 

B = £2si-nh^il3T), (29) 
C = 4sinh^(/3T), (30) 
where €i = X , £2 = £o£i, £3 = and K is an integrating constant. 

Therefore, the metric ^ reduces to 

ds^ = -dt^ + il sinh5^ {(3T)dx^ + £\ sinh^ {PT)e-^°"'dy'^ 

+ ^3sinh^(/3r)dz2. (31) 

In the following sections we deal with two cases, (i) non-interacting two-fluid model and (ii) interacting two- 
fluid model. 

4 Non-Interacting Two-Fluid Model 

In this section we assume that two-fluid do not interact with each other. Therefor, the general form of conser- 
vation equation (jl7p leads us to write the conservation equation for the barotropic and dark fluid separately 
as, 

+ 3" fpM + pi^)) = pM + (1 + J"'))p("'\2n + 1)§ = 0, (32) 

and 

p^"*^) + 3- + p'^'^'A = p'^-^"^ + (1 + J-^''^)p^'^''\2n + 1)5 = 0. (33) 

a \ J 6 



Integration of ((32l) leads to 

p(-) = ppC-(2"+i)(i+"""') = po^osinh''^(i+"""')(/3T), (34) 
where po is an integrating constant and lo = £^ (2n+i)(i+w ) ^ 
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Figure 1: The EoS DE parameter oj^'*'^) versus T for 
n = l3 = a = £3 = Iq = 1, po = 10. 



Figure 2: The total energy density parameter fl vs. T 

for n = = a = £3 = 1 



By using Eqs. (HH), ([201) and dM]) in Eqs. (US]) and ([12]), we obtain 



+ polo sinh-3(i+"""') (^^^^ 



and 



2n- + n(3n-2)— - — 



w(™Vo^osinh-3(^+'"""')(/3T). 



(35) 
(36) 



Using Eq. ^ in Eqs. ([35]) and (|36l), we obtain the values of p'-'^''^ and p^'^'^^ as 

= ^^if coth^(/3T) - a^^3-2"sinh-T^(/3r) - po^o sinh-3(i+"""') (/3T), (37) 

oj(™Vo^osinh-3(i+'"""''(/3T). (38) 



^("'+^^+^)coth^(/3T)-%±l)^co.ec/.^(/3r) 



(2n + 1)2 



(2n+ 1) 



respectively. 

The EoS (oj'^'^'^') of DE, for model pip is found to be dark energy in term of time as 



gggffi) coth^(/3T) - 3M^eo.ec/z^(/?r) +c.(")po^osinh'^(^-""')(/3r) 
"^^;;+gf coth2(/3T) - a2f^2ng.^j^-(^ (^y) _ pp/„sinh-3(i+"<'"')(y3T) 



(39) 



The behaviour of EoS for dark energy in term of cosmic time T is shown in Fig. 1. We observe from the Fig. 
1 that for w^™) > 1, uj^'^'^'> varies from non-dark region crossing quintessence region and ultimately approaching 
to cosmological constant region (w*^'*'^) = —1). But for o;'^™) < 1, the variation of w^'*'^) starts from from super 
phantom region (w™ > 1, w^"*^) < -3) crossing phantom region (w^'*^) < ~1) and finally approaches to cosmolog- 
ical constant region (cj^'*'^-' = — 1). The EoS parameter of the DE may begin in non-dark region or quintessence 
{uj"^" > -1) region and tends to —1 (cosmological constant, uj'^'^ ~ —1) by exhibiting various patterns as T 
increases; see Fig. 1. On the other hand, while the current cosmological data from SNIa (Supernova Legacy 
Survey, Gold Sample of Hubble Space Telescope) ^M\- CMB (WMAP, BOOMERANG) [Sll |97] and large 
scale structure (SDSS) [68j data rule out that w^"*^) ^ — 1, they mildly favour dynamically evolving DE crossing 
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the PDL (see [98—100, 57]) for theoretical and observational status of crossing the PDL). Thus, our DE model 
is in good agreement with well established theoretical result as well as the recent observations. 




and 



Jd. 



Figure 3: The plot of energy conditions versus T for Figure 4: The plot of deceleration parameter q Vs. T 
n = /3 a = ^3 = /o = 1, Po = 10, ujm = 0.5 for a = 1 

The expressions for the matter-energy density fl^"^^ and dark-energy density 57^'^'^) are given by 

^(™) ^P^^ po/osinh-^(^-""')(/?T) 

3if2 3/32 coth^ (/3T) ' ^ ) 

u,, P-^) ^fZV^f coth^(/3T) - a^^3-^"sinh-™(/3r) - po/o sinh-^(^+"""') (/3T) 

Vl^ > = TT = — n . (41) 

3iJ2 3/32 coth^ (^T) ^ ' 

respectively. Adding Eqs. (^(1)) and (pij) . we obtain total energy (fi) 

, , M!i+H)|!coth2(/3r)-a2£-2«sinh-2lTT(/3T) 

3/32 coth^ (/3r) ^ ' 

Figure 2 depicts the variation of the density parameter {VL) versus cosmic time T during the evolution of the 
universe. From the Figure 2, it can be seen that the total energy density f2 tends to 1 for sufficiently large time 
which is supported by the current observations. 

The dark energy with w''''^' < —1, the phantom component of the universe, leads to uncommon cosmological 
scenarios as it was pointed out by Caldwell et al. [101] . First of all, there is a violation of the dominant energy 
condition (DEC), since p^'^"'' + p<-d-^) < 0. The energy density grows up to infinity in a finite time, which leads 
to a big rip, characterized by a scale factor blowing up in this finite time. These sudden future singularities 
are, nevertheless, not necessarily produced by a fluid violating DEC. Cosmological solutions for phantom matter 
which violates the weak energy condition were found by Dabrowski et al. |102| . Caldwell [103) . Srivastava [24] . 
Yadav [104] have investigated phantom models with w^''"^) < — 1 and also suggested that at late time, phantom 
energy has appeared as a potential DE candidate which violets the weak as well as strong energy condition. 

The left hand side of energy conditions have been depicted in Figure 3 for different values of T. From Figure 
3, for w'™-' < 0.5 (i.e. phantom model) (also see Figure 1) , we observe that 

(i) p'-^"^ > 0, (m) p^'^^^ + p'-'^"^ < 0, {tit) p^'^'^ + ip'^^"^ < 0. 
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Thus, from above expressions, we observe the phantom model violates both the strong and weak energy condi- 
tions, as expected. 



The expressions for physical parameters such as directional Hubble parameters {H^, Hy, H^), the Hubble 
parameter (H), scalar of expansion (0), shear scalar (cr), spatial volume V, the anisotropy parameter and 
the deceleration parameter [q) are, respectively, given by 



3n 



2n+ 1 



coth {/3T), 



2n + l 



coth(/3r). 



= 3H ^ 3/3 coth {I3T), 

2 



1 



1 



/32 coth^ (^T), 



V = sin^ (/3T)e" 



1 



(43) 

(44) 
(45) 
(46) 
(47) 
(48) 



q = -tanh^ (/JT). 



(49) 

It is observed that at T = 0, the spatial volume vanishes and other parameters 9, a, H diverge. Hence the 
model starts with a big bang singularity at T = 0. This is a Point Type singularity [105] since directional scale 
factor A{t), B{t) and C{t) vanish at initial time. For n = 1, Am = 0. Hence the model is isotropic when n = 1. 
Fig. 4 shows the evolution trajectory of deceleration parameter, from which we observe that the evolution of 
the universe is in accelerating phase at present era. This is good agreement with recent observations that our 
universe is in an accelerated expansion moment [T] — [3]- 

We also solve for the deceleration parameter g as a function of the redshift z = — 1 + where Oq is the 
present value of the scale factor (i.e at z = 0). It is given by 



tanh < sinh 



-90 



(go + l)(2 + l)^ 



(50) 



Here go is the present value of the deceleration parameter i.e. at z = 0. The analyses of cosmological observations 
in literature [106) — [TTT| furnish the transition redshift of the accelerating expansion as given by 0.3 < zt < 0.8. 



5 Cosmic Jerk Parameter 

A convenient method to describe models close to A CDM is based on the cosmic jerk parameter j, a dimensionless 
third derivative of the scale factor with respect to the cosmic time |112) — [TTB] . A deceleration-to- acceleration 
transition occurs for models with a positive value of jo and negative qq. Flat A CDM models have a constant 
jerk j = 1. The jerk parameter in cosmology is defined as the dimensionless third derivative of the scale factor 
with respect to cosmic time 

Mt) = (51) 

where the 'dots' denote derivatives with respect to cosmic time. The jerk parameter appears in the fourth term 
of a Taylor expansion of the scale factor around ao 

^ = l + Ho{t- to) - \qoHl{t - tof + poH^t - to)' + 0[{t- to)'] , (52) 
ao ZD 

where the subscript shows the present value. One can rewrite Eq. (jSlT) as 

j{t) = q + 2q^-j^. (53) 
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Using Eq. ^SU^) in Eq. ([S3]) we obtain 



Jit) 



3(cosh^ (/3r) - 1) 



(54) 



cosh^ (/3T) 

This value overlaps with the value j ~ 2.16 obtained from the combination of three kinematical data sets: the 
gold sample of type la supernovae [HS], the SNIa data from the SNLS project [55], and the X-ray galaxy cluster 
distance measurements [117] for 

cosh {VS) 



T 



f3 



(55) 
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Figure 5: The EoS DE parameter oj^'^'^) versus T for Figure 6: The EoS DE parameter uj^'^'^^ versus T for 
n = /3 = a = 4 = I = 1, po = 10, K = 0.03 n = ^ = a = 4 = / = 1, Po = 10, = 0.5, k = 0.03 



6 Interacting Two-Fluid Model 

In this section we consider the interaction between dark and barotropic fluids. For this purpose we can write 
the continuity equations for barotropic and dark fluids as 



p(")+3- fp(™)+p(™)) =p(™) + (l+c^('"))p('")(2^+l)S =Q 
a V / C 



and 



^(de) 



^-0 

a V 



{de) 



(1 



(de) 



)p(''^)(27' 



(56) 



(57) 



The quantity Q expresses the interaction between the dark components. Since we are interested in an energy 
transfer from the dark energy to dark matter, we consider Q > 0. Q > 0, ensures that the second law of 
thermodynamics is fulfilled |118] . Here we emphasize that the continuity Eqs. ([5S)) and ([57)1 imply that the 
interaction term (Q) should be proportional to a quantity with units of inverse of time i.e Q oc j. Therefor, 
a first and natural candidate can be the Hubble factor H multiplied with the energy density. Generally, Q 
could be taken in any of the following forms: (i) Q cx 1I20], (ii) Q cx Hp^'^'='>^^ [m], or (iii) 



P 



(m) 



) |122j - jl24j . It is worth to mention that because of our lack of knowledge of the nature of 



DE and DM, one can choose any of the above specific form of the interaction term Q. Moreover, as noted in ref 
[85] ■ describing the interaction between the dark components of the universe through a microphysical model is 
not available today. In our study. Following Amendola et al. [119j and Gou et al. [120) . we consider 

Q = Hkp^"'\ (58) 
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where k is a coupling coefficient which can be considered as a constant or variable parameter of redshif. As 
it is shown by Guo et al. |115j . the combination of the SNLS, CMB, and BAO databases marginalized over a 
present dark energy density gives stringent constraints on the coupling, —0.08 < k < 0.03 (95% C.L.) in the con- 
stant coupling model and —0.4 < kq < 0.1 (95% C.L.) in the varying coupling model, where kq is a present value. 



Using Eq. ([55)1 in Eq. ([55]) and after integrating, we obtain 



po/sinh-3(i+'^*""-'=) (^y^^ 



(59) 



where I — (2n+i)(i+'^ ^q^q ihdX here for simplicity a constant coupling constant, k, is considered. 

By using Eqs. (HH), (gHl) and ^ in Eqs. ^ and we obtain 

^^^-(2n+l)(l+a;<'"'-K)^ 



and 



pile) 


= n{n 




C<2n 




h n{3n 
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C2n 



PO 



Using Eq. p6| in Eqs. (|60l) and (iGlT) . we obtain the values of p^''^) and as 



^ 9n(n + 2)^ coth^(/3T) - a^fg-^" sinh'^^ (/3T) - polsmh-^^'+"'"^'-'^'> (f^T), 
(2n + 1)^ 



(60) 
(61) 

(62) 



and 



6n^2 3n(5n-2);92 



(2n+l) (2n+l)2 

-po(^^"'^ -K)^sinh-3(i+"' 



coth^ {I3T) - a2^;^2"sinh^^ (^T) 



(63) 



respectively. 

Also the EoS parameter for DE (w*^'*'^)) is obtained as 



(de) 



(fgfi) + 'tn+if (/^^) - «^^3~'"sinh-^ (/3T) +po(^('") - ^)^sinh-^(^+"""'-) 

""lilX^^f coth2(/3T) - a2^-2nsini,-^ (^7^) _ pn/sinh-3(i+"<'"'-«) ^^^^ 



(/3r) 



(64) 

The behavior of EoS (w'-'^'^') for dark energy in term of cosmic time T is shown in Fig. 5. It is observed that 
the EoS parameter is an increasing function of time and the rapidity of its increase at the early stage depends 
on the value of w'-™-', while later on it tends to the same constant value (i.e, —1) independent to it. 

The expressions for the matter-energy density ri'^™^ dark-energy density Vl'^'^'^^ and the density parameter f2 
are same as given by Eqs. (HUl) - in the non-interacting case. 

Based on Eqs. (15^ and (1551) . the left hand side of energy conditions have been plotted in Figure 6 for different 
values of T. From Figure 6, for w^'"^ < 1 (i.e. phantom model) (also see Figure 5) , we observe that 

(«)/5^'^'^>0, (m) p^'^'^) -Hp^"*") < 0, (zm) -f 3p('*'=) < 0. 

Thus, from above expressions, we observe the phantom model violates both the strong and weak energy condi- 
tions, as expected. 

Studying the interaction between the dark energy and ordinary matter will open a possibility of detecting 
the dark energy. It should be pointed out that evidence was recently provided by the Abell Cluster A586 in 
support of the interaction between dark energy and dark matter |125[ 1126) . We observe that in non- interacting 
case both open and flat universes can cross the phantom region whereas in interacting case only open universe 
can cross phantom region. 
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Figure 7: The variation of s against r 



Sahni et al. |127j . Alam et al. [27] have introduced a pair of parameters {r, s}, called Statefinder parameters. 
In fact, trajectories in the {r, s} plane corresponding to different cosmological models demonstrate qualitatively 
different behaviour. The Statefinder parameters can effectively differentiate between different form of dark 
energy and provide simple diagnosis regarding whether a particular model fits into the basic observational data. 
The above Statefinder diagnostic pair has the following form: 



H H , r-1 
For our model, the parameters {r, s} can be explicitly written in terms of T as 



2 1 
tanh (aT) , s — - 
o 



sech'^ (aT) 
tanh^(aT) + i 



So the relation between r and s has the implicit form; 

3(2r + l)s + 2(r- 1) 0. 



(65) 



(66) 



(67) 



From Fig. 7, we observe that s is negative when r > 1. The figure shows that the universe starts from an 
asymptotic Einstein static era (r oo, s — oo) and goes to the ACDM model (r = 1, s = 0). 



7 Concluding Remarks 

In this paper, we have studied a spatially homogeneous and anisotropic Bianchi type-Ill space time filled with 
barotropic fiuid and dark energy possessing dynamic energy density. The role of two-fluid cither minimally or 
directly coupled in the evolution of the dark energy parameter has been investigated by considering the scalar 
expansion in the model is proportional to the shear scalar. The held equations have been solved exactly with 
suitable physical assumptions. To prevail the deterministic solution we choose the scale factor a{t) = sinh(/3r), 
which yields a time-dependent deceleration parameter, representing an accelerating model. The solutions satisfy 
the energy conservation Eq. pT|) and the Raychaudhuri Eq. (|18l) identically. Therefore, exact and physically 
viable Bianchi type-Ill model has been obtained. The main features of the model are as follows: 

• For n = 1 the anisotropic parameter Am tends to zero. Hence, the present model is isotropic at n = 1. 
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• The derived DE model represents an acceleration universe (see, Figure 4) which is in good agreement with 
recent observations [J — [3]. 

• In non-interacting two- fluid model, we observe that for w^™) > 1, uj^'^"') varies from non-dark region cross- 
ing to quintessence region and ultimately approaches to cosmological constant region 

(^(de) = But for 

< 1^ the variation of o;'^'*'^^ starts from super phantom region crossing phantom region 

(t^(rfe) < _i) and 

finally approaches to cosmological constant region 

• In interacting two-fluid model, we observe that uj^'^'^'^ is rapidly increasing function of time and ultimately 
approaches to cosmological constant region 

• In non-interacting and interacting two-fluid models, we obtained that phantom model violates both the 
strong and weak energy conditions, as expected (see, Figures 3 & 6). 

• In both non-interacting and interacting two-fluid scenario, the total density parameter (Q) approaches to 
1 for sufficiently large time (see. Figure 2) which is reproducible with current observations. 

• The cosmic jerk parameter in our descended model is also found to be in good agreement with the recent 
data of astrophysical observations namely the gold sample of type la supernovae [96j, the SNIa data from the 
SNLS project [9^, and the X-ray galaxy cluster distance measurements |118j . 

• In both non-interacting and interacting two-fluid scenario, it is observed that such DE models are also in 
good harmony with current observations. Thus, the solutions demonstrated in this paper may be useful for bet- 
ter understanding of the characteristic of anisotropic DE in the evolution of the universe within the framework 
of Bianchi type-Ill space-time. 

• The Stateflnder pair {r, s} enable the behaviour of different stages of the evolution of the universe i.e. the 
universe starts from asymptotic Einstein static era (r oo, s — >■ — oo) and goes to A CDM model (r — 1, s — 0). 

Finally, the solutions presented in this work can be one of the potential candidates to describe the observed 
universe. 
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